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Abstract 


In this paper the general problem of a tapered specimen containing an 
edge crack is formulated in terms of a system of singular integral equations. . 
The equations are solved and the stress intensity factor is calculated for 
a "compact" and for a "slender" tapered specimen, the latter simulating 
the double cantilever beam. The results are obtained primarily for a pair 
of concentrated forces and for crack surface wedge forces. The stress 
intensity factors are also obtained for a long strip under uniform tension 
which contains inclined edge cracks. 


1. Introduction 

The main objective of this paper is to provide a solution for a tapered 
specimen containing an edge crack by using an integral equation technique. 

The geometry of the specimen is shown in Fig,’ 1 where the dimensions H, B, 

£, m, n, and e are arbitrary. By choosing the relative dimensions properly 
one may simulate either the tapered compact tension specimen (Fig. 1) or 
the tapered double cantilever beam specimen (Fig. 2 ). Even though the results 
are given only for two types of loading shown in Figures (la) and (lb), the 
formulation of the problem is quite general and any arbitrary loading on 
the crack surfaces or on the inclined boundaries, for example, can easily 
be taken into account. Of course, the problem can also be solved by using 
a finite element method. Because of its practical importance the rectangu- 
lar specimen having an edge crack has been studied rather widely (see, for 
example, [1] and the articles by Bowie, Isida, and Wilson in [2]), The 
general problem was also discussed in a recent article [3], 
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2. Formulation of the Problem 

We will first consider the general crack problem for an infinite strip 
described in Fig. 3. It will be assumed that y=0 is a plane of symmetry 
and the location (nun) of the concentrated force P and the crack surface 
tractions are arbitrary. For the purpose of deriving the integral equations 
one may express the stress state at a point (x,y) in the strip as follows: 

a^Cx.y) = a-iijU.y) + CT 2i + CT 3ij( x ’ y ) + G 4ij( x ’ y ) 

+ a 5i j(x,y) , (i »j=x,y) . (1) 


where the stress components are associated with an infinite plane con- 
taining a crack along (y=0, a<x<b), a 2 .jj and are associated with a plane 
having a crack along (c<r<d, 0 ) and (c<r<d, - 0 ), respectively, o^- ■ relate 

to an infinite plane without any cracks under concentrated forces P at x=m, 
y=+n, and a^. . are associated with an infinite strip. The total stress 
state cr.j j must satisfy the following boundary conditions (Fig. 3): 


<T x) .(0,y) 

II 

O 

q 

3 

(Q,y) = 0, 

(0<y<») , 

(2) 

CT xx (W - y) 

J? 

•t> 

A 

O 

II 

(H,y) = 0, 

( 0 <y <«» ) , 

(3) 

°xy (x,0) 

= 0, (0<x<H) , 


(4) 

Oyytx.°) 

n 

"x 

CL 

II 

(a<x<b) , 


(5) 

°ee (r ’ e > 

= P 2 (r), 

(c<r<d) , 


(6) 

o r 0 ( r > 6 > 

■ P 3 (r) > 

(c<r<d) ' , 


(7) 


where p^ , p 2 and p 3 are known functions and where it is assumed that because 
of symmetry only one half (y>0) of the domain needs to be considered. 

We now define the following unknown functions: 

g-jOO = [Uy(x,+0)-Uy(x,-0)], (a<x<b) , (8) 


9 2 (r) = C u 0 (r,9+O)-u 0 (r,e-O)], (c<r<d) , 


(9) 
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g 3 (r) = [u r (r,e+0)-u r (r,6-0)] f (c<r<d) , 



(10) 

g 4 (r) = Cu 0 (r,-e+O)-u 0 (r,-e-O)], (c<r<d) , 



(ID 

g 5 (r) = jjp- [u r (r,-e+0)-u r (r,-e-0)], (c<r<d) 



(12) 

where u , u or u , u are the components of the displacement 

A Jr I U 

to x,y or r,e coordinates, respectively (Fig. 3). The stress 
ated with an infinite plane containing a symmetrically loaded 
x axis may then be expressed as 

k 

vector referred 
state associ- 
crack on the 

u 

a li j(x,y) = | G. j .(x,y,t,0)g 1 (t)dt , (i,j=x,y) 
a 



(13) 


where the kernels G. . are given in Appendix A. The stresses associated 
with edge dislocations of densities g 2 and 93 distributed along (c<r<d, e = 
constant) are given by 


d 


c 2ij (x,y) = j 

u 

G . j ( x ,'y ,x Q ,y Q ) [g 2 ( r Q ) cos e+g 3 ( r Q ) si n e] dr Q 
c 


- 

d 

• 


+ 

• 

H ij^ x,y ’ x o ,y o^ 9 3^ r o^ cose “ 9 2^ r o^ sine ^ dr o ’ 
c 



(i »j=x,y) , 

(14) 

where G i ^(x.y ,x Q ,y o ) and H i . (x^,x Q ,y 0 ) , (i,j=x,y) are again given in 
Appendix A, and 


x_ = r cose . 
0 0 

. y 0 = B + i^sine . 

(15) 


Similarly, for the crack lying along (c<r<d, - 0 ) we have 

d 

CT 3i j ( x,y ) = { G ij^ x,y,x o’ y o^ 9 4^ r o^ cos6 " 9 5^ r o^ sin0 -' dr o 


+ j H. J .(x,y,x o ,y o )Cg 5 (r o )cos 0 +g 4 (r o )sine]dr o , (i,j=x,y), 
c 

- 3 ~ 


(16) 



(17) 


x 0 = r 0 cose , y Q = -B-r Q sine . 


The stresses due to a pair of concentrated forces P and -P (per unit 
thickness) are given by (Fig. 3) 


a 4ij -(x,y) = P Qj .(x,y,m,n) , (i,j=x,y) 


08) 


where the functions Q. . are given in Appendix B. 

• J 

Finally, by using Fourier transforms and the symmetry of the problem, 
the stresses in an infinite strip 0<x<H, -®<y«» may be expressed as 


00 

a 5xx (x,y) = - j- {[ a (A 1 +xA 2 ) +~- Agje' 


aX 


1+tc 


+ Ea(A 3 +xA 4 ) --y 1 A 4 Je aX }cosayda 


a 5yy (x,y) = IT 


{[a(Ai+xA 2 ) + A 2 ]e 


■aX 


k-3 


+ [a(A 3 +xA 4 ) -•^ 2 ^ A 4 Je aX }cosayda 


W x ’ y) = ' T 


K-l 


{[a(A^+xA 2 ) + -~ 2 ~ A 2 ]e 


-ax 


<-l 


- [a(A 3 +xA 4 ) - A 4 ]e aX }sinctyda 


(19a-c) 


where A^,...,A 4 are unknown functions of a. 

If we now substitute from (13)-(19) into (1), use the homogeneous boun- 
dary conditions (2) and (3), invert the Fourier transforms and evaluate the 
resulting infinite integrals, we obtain the unknown functions A-|,...,A 4 in 
the following form: 


A- (a) = l 
1 k=l 


B ik (a»' t ) g k( t ) dt + PC i(“) » (i =1 » — >4, 


a-j=a, b^=b, a k =c, b k =d, k=2,3) , 


( 20 ) 


- 4 - 



where the symmetry properties 


u Q (r,-e+O)-u 0 (r,-e-O) = u 0 (r,9+O)-u 0 (r,e-O), 


u r (r,-e+0)-u r {r,-e-0) = -[u r (r,e+0)-u r (r,e-0)] (21a,b) 

have been used. The explicit expressions of A* may be found in Appendix C. 

We now observe that the boundary condition (4) is satisfied by the 
assumed symmetry and from (9 )- (12 ) and (21) we have 


9 4 (r) = g 2 (r) , g 5 (r) = -g 3 (r) 


( 22 ) 


The three remaining boundary conditions (5)- (7) may then be used to deter- 
mine the unknown functions g-j, g 2 anc * 93 * Thus, by substituting from (13)- 
(20) into (5 )- (7) and by using (22) we obtain the following system of inte- 
gral equations to determine g-j , g 2 and g 3 


3 

E 


j=l 


b J 


a j 


h^ j(s,t)gj (t)dt = p..(s) + P h. (s) , (i=l ,2,3; 
a-|=a<s=x<b=bi , a^=c<s=r<d=b| < , k=l,2) , 


(23) 


where h.., (i ,3=1 ,2,3) and h. , (i =1 ,2,3) are known functions and may be 
expressed in terms of infinite integrals of the following form by using 
the information given in the appendices A, B and C: 


h ij -(s,t) = 


K..(s,t,a)da , (i ,j=l ,2,3) 


h i (s) 


Kj (s ,a)da , (i=l ,2,3) . 

0 


(24) 

(25) 


For the crack problem shown in Fig. 3, by separating singular part of the 
kernels h. . through an asymptotic analysis, it can be shown that the main 

* vJ 

diagonal elements of h. . have Cauchy type singularities. That is, for 

* 0 

a-H*> if we let 
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( 26 ) 


K.j j (s ,t,a) -*■ K ^ j oo ( s , 1 ,<x ) 
we find 


K- - ( s ,t,a)da 
j 11“ 

o 


2 V ] 

7T ( 1 +K ) t“S * 


(i=l ,2,3) 


(27) 


and the kernels h. . may be expressed as 

* J 


h ij (s - t) = 


2v 


5.. 


tt(1+k) 


tlj+ k. .(s,t) , (i ,j=l ,2 ,3) 


U 


(28) 


where the functions are bounded within the closed interval aj_<(s ,t)<b^ , 
(j=l,2,3). From the definitions given by (8)- (10 ) and from Fig. 3 it is 
clear that the density functions g.(t), (j=l,2,3) must satisfy 

J 


b 

g,(t)dt = 0 , 
J 1 
a 


d 

g k (t)dt = 0 

J 

c 


(k=2,3) 


(29) 


Referring to, for example, [4], it is known that the solution of the 
system of singular integral equations may be expressed as 


g^t) = 
= 


F,(t) 


/(t-a)(b-t) 
/(t-c) (d-t) 


, (a<t<b) , 

, (k=l,2; c<t<d) 


(30) 


(31) 


The Modes I and II crack tip stress intensity factors may then be defined 
by and evaluated from the following relations: 


k-j (a) = 

= lim 

/2(a-x) 

a yy (x,0) 

_ 2y 
1+K 

lim /2(x-a) 

g-j(x) , 




x-va 


x-»a 




k-| (b) = 

= lim 
x->b 

-Q 

1 

X 

ffyyU.O) 

= - lim /2(b-x 
1 K x-*b 

7 9] (x) , 

(32 

a,b) 

k-j ( c) = 

= lim 

i\3 

n 

i 

-j 

°ee (r ’ e) 

= iL. 

1+K 

lim /2(r-c) 

g 2 ( r ) » 



r-*c 


r->c 



k 2 (c) = 

s lim 
r-*c 

/2(c-r) 

o re (r .e) 

= Jh_ 

1+K 

lim /2(r-c) 
r-K: 

g 3 (r) , 

(33 

a,b) 
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(34 a,b) 


k-j(d) = lira /2(r-d) a Q0 (r,e) = - ^Mim /2(d-r) g 2 (r) , 

r-d r-»d 

k 2 (d) = Urn /2(r-d) a (r,e) = - lim V2 (d-r) g 3 (r) . 

r-*d r-vd 

The functions h.(s), (1=1 ,2,3) and the kernels h..(s,t), (i,j=l,2,3) 

I I J 

may be evaluated in terms of G. ., H. ., Q. . , (i,j=x,y), and A- , (i=l,...,4) 

» J • J I J t 

in a straightforward manner. However, the manipulations and the resulting 
expressions are quite lengthy and will not be reproduced in this paper. 

We also note that in obtaining the integral equations from the crack 
surface boundary conditions (6) and (7), the stress state cr--(x,y), (i ,j= 

* J 

x,y) as given by (1) is substituted into 

P 2^ = a 9e^ r ’ e ) = ^XX (X ’ y)Sin20 + ^yy( X »y) COs2e - 2t7 X y( X >y) S ' in 0 COS0 > 
P 3 ( r ) = a r9 (r,e)=[ff yy (x,y)-a xx (x,y)]sinecose +a xy (x,y)(cos 2 e-sin 2 e) , 

x = rcose , y=B+rsine , (e=constant, c<r<d) . (35 a-d) 

For the case of "edge" cracks, i.e,, for a=0 or c=0, the asymptotic 
analysis would show that the singular parts of the kernels h..(s,t) are 

* J 

generalized Cauchy kernels. That is, as shown, for example, in [5] the 
singular parts contain, in addition to the standard Cauchy kernels, terms 
that become unbounded as the variables s and t approach zero simultaneously 
and as a consequence at the end point t=0 the density functions g.(t), 
(i=l,2,3) become bounded (see, also [3]). The integral equations were 
solved by using the technique described, for example, in [6] (see also, 

[5]). 


3. Numerical Results 

The first numerical example considered is described by the insert in 
Fig. 4, namely a long strip containing two symmetrically located inclined 
edge cracks. Rather extensive results for parallel cracks may be found in 
[3]. For a given B/H ratio (of 0.2) the normalized stress intensity factors 
are given in Table 1, In this example, the relative crack length Z and the 
angle 9 are the variables. Table 1 shows the results for a uniform tension 
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°yy =a o away ^ ronl the cractc region. In this case the crack surface trac- 
tions for the perturbation problem are the input functions in the integral 
equations (23) and are given by 

P 2 (r) = -cr o cos 2 e, P 3 (r) = -c o sin8cose . (36) 

The special case of this problem for 8=0 was studied in [3] and are also 
calculated here for comparison. It may be observed that for small values 
of e k 2 (d) is negative and, as pointed out in [33, would force the cracks 
to grow away from each other. However, for relatively large values of e 
k 2 (d) becomes positive. This means that for such angles the cracks tend 
to grow toward each other or, they tend to orient themselves more nearly 
perpendicular to the direction of the external load. This, of course, is 
the physically expected result. Theoretically, if the Mode II stress inten- 
sity factor k 2 (d) is zero, momentarily the crack would be expected to grow 
in its current plane. Table 1 implies that for a given crack length 1 / H 
the value of e=9 0 corresponding to k 2 (d)=0 can be calculated. For B/H = 

0.2 the calculated values of e„ are shown in Fig. 4. 

Table 2 shows the results for three edge cracks of equal length in a 
long strip under uniform tension a 0 away from the crack region. In this 
case too the results for 8=0 agree with those given in [3], 

Tables 3-5 show the results for a tapered "compact" specimen contain- 
ing an edge crack (a=0, b=£), The stress intensity factor given in Table 
3 correspond to a uniformly pressurized crack shown in Fig, (la) for B/H = 
0.2 and B/H = 0.48, Tables 4 and 5 give the stress intensity factors for a 
pair of concentrated forces and crack surface wedge forces, respectively. 

The stress intensity factor for a "slender" tapered specimen (or the 
tapered double cantilever beam specimen) described in Fig. 2 is given in 
Fig. 5. The figure indicates that for B/H=0.1 one could have a "constant 
k" regime only for 8=30° and 0.15<£/H<0.4. For the double cantilever beam 
specimen an approximate value of the stress intensity factor may also be 
obtained by using the following energy balance relation (Fig. 2): 

a -rfk^b) 

■k (u - v) * -T- • (37) 
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where L is the length of the "beam" (Fig. 2), U is the work of the external 
forces P, V is the strain energy, E'=E for plane stress, and E'=E/(l-v 2 * ) 
for plane strain. By using the beam theory and by taking into considera- 
tion the energy due to the transverse shear as well as the bending stresses, 
for the specimen shown in Fig. 2 the stress intensity factor may be esti- 
mated as follows: 

Mb) = /F (38) 

0 

where the shape factor is calculated to be 

f . J=l + Hi 

*7 5 

for 9=0 and 

h , , h 2 (h +2Ltane) 

r _ 0 / 1+V , 1 \ _ 0__0 

h Q +Ltane V ~1T~ tan 7 ? 7 " tan^e^+Ltane)* 

for 8>0 (Fig. 2). A limited comparison of the stress intensity factors 
obtained from the elasticity and the beam theories is shown in Table 6. 
Considering the simplicity of the beam solution, one may observe that the 
agreement between the two results is fairly good. However, one may also 
observe that the beam theory gives consistently smaller values than that 
obtained from the elasticity solution. Physically, this is not really sur- 
prising. As Fig. 2 indicates, the constraint of the actual specimen is 
less than that of the pair of cantilever beams. Consequently an end con- 
dition which is more realistic than the built-in end of a cantilever beam 
would give a greater stress intensity factor. 
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Table 1. Modes I and II normalized stress intensity factors in a long strip 
under uniform tension cr q which contains two syirnietrically located 
inclined edge cracks, B/H = 0.2, c=0, a=b (Fig. 3). 



e 

m 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 



BBI 

IB 

IB 

IB 

2.5980 

3.8235 

Md) 


UB 

IB 

IB 

IB 

2.6960 

3.9936 

1 

10° 

■ 

1.2358 

1.5128 

1.9647 

2.7002 

3.9974 

V* 

15° 

■I 

1.2169 

1.4870 

1.9210 

2.6201 

3.8559 


30° 

IB 

1.0521 

1.2531 

1.5591 

2.0314 



0° 


-0.1083 

-0.1694 

-0.2280 

-0.2937 

-0.3597 

k ? (d) 

'5° 

0.0319 

-0.0266 

-0.0642 

-0.0854 

-0.0880 

-0.0463 

L 

10° 

0.0953 

0.0551 

0.0416 

0.0526 

0.1087 

0.2302 

a Q /l 

15° 

0.1560 

0.1333 

0.1422 

0.1875 

0.2822 

0.4631 


30° 

0.3032 

0.3188 

0.3697 

0.4603 

0.6047 

0.8460 
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Table 2. Normalized stress intensity factors in a long strip under 

uniform tension a 0 which contains three edge cracks of equal 
length; B/H =0.2, c=0, a=0, d=b=£ (Fig. 3). 


e 

l 

k 1 (b)/or 0 vT 

k-j (d)/a Q /£ 

k 2 (d)/a 0 /e 


0.01 

1.1100 

1.1149 

0.0005 


0.11 

0.7923 

0.9849 

-0.1080 

0° 

0.21 

0.8311 

1.1134 

-0.1740 

0.31 

0.9929 

1.3656 

-0.2398 


0.41 

1.2902 

1.7747 

-0.3251 


0.51 

1 . 8205 

2.4480 

-0.4452 


0.01 

1.1102 

1.1091 

0.0592 


0.11 

0.8089 

0.9987 

-0.0485 

E 0 

0.21 

0.8618 

1.1394 

-0.1031 

0 

0.31 

1.0481 

1.4080 

-0.1462 


0.41 

1.3874 

1.8432 

-0.1950 


0.51 

1.9915 

2.5624 

-G.2485 

• 

0.01 

1.1112 

1.0625 

' 0.1727 


0.11 

0.8566 

0.9815 

0.0745 

15° 

0.21 

0.9466 

1.1265 

0.0388 

0.31 

1.1886 

1.3947 

0.0283 


0.41 

1.6102 

1.8251 

0.0355 


0.51 

2,3360 

2.5372 

0.0764 


0.01 

1.1136 

0.9145 

0.3049 


0.11 

0.9488 

0.8573 

0.2227 

co 

o 

o 

0.21 

1.0913 

0.9707 

0.2113 

0.31 

1.3933 

1.1781 

0.2393 


0.41 

1.8832 

1.5011 

0.3059 


0.51 

2.6784 

2.0131 

0.4375 
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Table 3. Normalized stress intensity factor lq(b)/cr 0 /£ in a 
tapered specimen containing an edge crack and sub- 
jected to uniform crack surface pressure (Fig, la). 



m 

B/H 

0 

0.1 

0.3 

0,5 

0.7 


0° 

IBB 

in 

3.0924 

6.4152 

0.48 

10° 

mm 

H 

2.9123 

6.3642 


mm 



2.8625 

6.3565 


H 

UB 

H 

2.8435 

6.3549 


0° 

mm 

4.0862 

7,3316 

10.7602 

0.2 

10° 

119 

2,7726 

3.9822 

.6.7785 


20° 

■ESI 

2.2798 

3.2527 

6.4318 


30° 

HI 

2.0052 

2.9977 

6.3694 
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Table 4. Normalized stress intensity factor k-|(b)/(P/E/H) 

in a tapered specimen which contains an edge crack and 
is subjected to a pair of concentrated forces P, a=0, 
B/H=0 .48, m/H=0.2, n/H=0.32, (Fig. lb). 



m 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0° 

2.8736 

3.6792 

4.5127 

5.2761 

6.2305 

7.8916 

sa 

2.8067 

3.4466 

4.1153 

4.8043 

5.8059 


E9 

2.7201 

3.3067 

3.9374 

4.6296 

5.6743 



Table 5. Normalized stress intensity factor ki (b)/(P/£/H) in 
a tapered specimen which contains an edge crack and 
is subjected to concentrated wedge forces P; a=0, 
B/H=0.48 , m/H=0.2, n=0 (Fig. lb). 



£7H 

C 

0.3 

0.4 

0.5 

0.6 

0.7 

0° 

5.3245 

5.4715 

6.2768 

7.8647 


10° 

4.9568 

5.0104 

5.8537 

7.5822 

11.2294 

20° 

4.8228 

4.8608 

5.7348 

7.5170 

11.2130 

CO 

o 

0 

4.7624 

4.7996 

5.6917 

7.4964 

11.2093 


Table 6. Comparison of the stress intensity factors k-, (b)/(P/>4T) 
calculated from the elasticity solution and from the 
beam theory for a slender tapered specimen; a=0, B/H= 
0.1, m/H=0.067, n/H=0.05, £-m=L, (Fig. 2). 


0 


L/H 


0.3 

0.4 

0.5 

0.6 

10° 

Elasticity 

9.70 

■61 


mm 

mwm 

Beam 

8.31 

■a 


Ks9 

Ha 

PO 

o 

o 

Elasticity 

6.38 

6.73 


6.89 

7.11 

Beam 

5.94 

6.37 


6.61 
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APPENDIX A 


The Green's functions G. . and H. . for a pair of edge dislocations 

* J * J 

at the point ( X Q , Y Q ) in an infinite plane. 
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APPENDIX B 


The Green's functions Q. . due to a pair of concentrated forces P 

* 0 

acting at the points (m,n) and (m,-n) in an infinite plane. 


Q xx (x ’ y ’ m ’ n) 2 tt( 1+k) *(x-m) 2 +(y-n} 2 " (x-m) z +(y-n) z ^ 


J^+n r.,_i 4(x-m) 2 i, (B1 ) 


r ■» H-y a-hi; 1 

( x-m ) z +(y+n ) z LK ~ " (x-m^+ly+n)^ 


Q yy (x,y,m,n) = 27 ^ { (x-ni^ ^(y-n )' 2 ' [ " (tc+3) + I3^F^y^F ] 


‘ rx-m)^+'(y+n ) 2 [ ‘ (lc+3) + Tx-S^y+n F ] } » 


(B2) 


Q xy (x,y,m,n) = ^(]+ K j- { ^mFT{ylnF [ ' (tc+3) + (x-m)*+(]-n) 2] 


" rx-ra'P+'(y+n) 2 ' I_(fC+3) + F-raP+Ty+n) 2 ^ * {B3) 
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APPENDIX C 


The functions A., (ct), (1=1 4). 

A- (a) = E E..(a)R.(ct) , (1=1,. ..,4) , (Cl) 

1 j=l 1J 3 

Ri (ct) = F.(a) + G.(a) + H.( a ) + P-(a) , (1=1,. ..,4) (C2) 

b 

F-|(a) = J C-j -j (a,t)g-j (t)dt , 

a 

d 

► 

G-j(a) = [C 12 (a,r 0 )g 2 (r 0 ) + C 12 (a,r 0 )taneg 3 (r 0 )]cosedr 0 , 

c 
.d 

H-j (ct) = CC 13 (cx,r 0 )g 3 (r 0 ) - C 13 (a,r 0 )taneg 2 (r o )]cos9dr ( ; , 

c 

P-j (ct) = C 14 (a)P , (C3) 

b 

F 2 (a) = Cg^ajtJg^tJdt , 
a 
, d 

S 2 (o) ■ C C 22*“’ r o* s 2^ r o* + C 22'“’ r o' tan09 3' I ’o^ COS9dr o ’ 

c 
d 

► 

H 2 (a) = [Cggfa.r* )g 3 (r ) " C 23^ a ’ r o^ taneg 2^ r o^ COSedr o ’ 

c 

P 2 («) = C 24 (a)P , (C4) 
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b 

F 3 (a) = C 31 (a.t) gi (t)dt , 
a 
f d 

63 (a) = C c 32^ a,r o^ 9 2^ r o^ + C 32^ ajr o^ tan9g 3^ r o^ COS0dr o » 

c 

,d 

H 3 (a) = [ c 33( a ’ r 0 )93( r 0 ) “ c 33 ^ a »»* o )tan0g 2 ( r Q )]cosedr 0 , 

c 

P 3 (a) = C 34 (a)P , (C5) 

b 

► 

F 4 (a) = C 41 (a,t)g-j (t)dt , 

a 

d 

G 4 (a) = [C 42 (a,r 0 )g 2 (r 0 ) + C 42 (a,r 0 )taneg 3 (r 0 )]cosedr 0 , 

c 

d 

H 4 (a) = [C 43 (cx ,r 0 )g 3 (r 0 ) - C 43 (a,r o )taneg 2 (r 0 )]cosedr 0 , 

c 

P 4 (a) = C 44 (a)P , (C6) 

C n (a,t) = - e' at , 

C 12^ a ’ r o^ = “ <+T ^ aX o e a ° cos ^o ) ’ 

c 13 (°* r o) = dr e " aX ° («V 1)sinay o * 

c 14 (a) = 4 e -am ( tc-1 -2am)si nan , (C7) 
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c 21 {a,t) = - e“ at (l-ot) , 

c 22 (a. r 0 ) = - e" aX ° (l-ax 0 )cosay 0 , 

C 23 (a ’ r o } = '^T aX o e " aX ° sinay o * 

c 24 (a) = j e" am (-<-l+2am)sinan , (C8) 

C 31 (o,t) = - ^ T c t (t-H)e‘ a(H " t) , 

C 32 (a ’ r o) = " «( x o -H)e“ a(H " x °^ cosay o , 

C 33^ a,r 0^ = i^T e " a ^ H " X0 ^a( H -X 0 )+ 1 ] sin ay 0 , 

C 34 (a) = 4y ~ ( ic+ ' l ) e ~ a(H " m) K-l-2a(H-m) sinan , (C9) 

C 41 («.t) = - e"“ (H_t) [l-a(H-t)J , 

C 42 (“ * r 0 ) = ” e " a(H " Xo) n-cx(H-x 0 )]cosay 0 , 

C 43 (a ’ r o ) = “dr «(x 0 -H)e- a(H " x o) sin a y 0 , 

C 44 (a) = e m ^[ic+l-2a(H-m)]sinan , (CIO) 

x o = r Q cos0 , y Q = B + r o sin 0 (Cll) 
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E-j i ( cx ) = (2ctD )[4 cc 2 H 2 -(k- 1 ) (-e 2 a ^-2ctH+l )] , 

Ei 2 (a) = (2ctD ) [4a 2 H 2 -(ic+"I ) (e 2 a ^-2aH-1 )3 > 

E 13 (a) = (2aD“ 1 )r(l-<)(e aH -e" aH ) + 2 a H(e' aH - K e aH )J , 

E u (a) = (2aD“ 1 )[(l+<)(e aH -e' aH )-2 a H(e" aH +,ce aH )] , (C12) 

E 21 (a) = D 1 (-e^ a ^-2aH+l ) , 

E 22 (a) = D' 1 (e 2 aH -2aH-l) , 

E 23 (a) = D ” 1 (2aHe aH +e ctH -e”“ H ) , 

E 24 (a) = D“ 1 (2aHe aH -e oH 4€“ aH ) , (C13) 

E 31 (a-) = (2aD 1 )[4c£ 2 H 2 +(k:- 1 ) (e 2 a ^-2aH-l)] , 

£ 32 ( 0 ) = (2aD ^)[(< + l)(e ^ a ^+2aH-1 )-4a 2 H 2 J , 

E 33 (a) = (2aD * )[(<-! )(e a ^-e a ^)+2aH(i<e a ^-e a ^)] , 

E 34 (a) = (2aD" 1 )C(<+l)(e aH -e’ aH )-2 a H(Ke“ aH +e aH )] , (C14) 

E 41 (a) = D~ 1 (e“ 2 “ H -2oH-l) , 

E 42 ( a ) = D" 1 (e" 2 ctH +2aH-l ) , 

E 43 (a) = D‘ 1 (e aH -e‘ aH +2aHe‘ aH ) , 

E 44 (a) = D" 1 (e a ^-e” aH -2aHe" a ^) , (C15) 

D(a) = e 2 aH +e' 2 aH -4a 2 H 2 -2 (C15) 
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(a) (b) 

Fig. 1 Geometry and loading conditions for tapered "compact" specimen 
with an edge crack. 
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